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Abstract 

Boltzmann equation describes the time development of the velocity 
distribution in the continuum fluid matter. We formulate the equation 
using the field theory where the velocity- field plays the central role. 
The matter (constituent particles) fields appear as the density and 
the viscosity. Fluctuation is examined, and is clearly discriminated 
from the quantum effect. The time variable is emergently introduced 
through the computational process step. The collision term, for the 
(velocity)**4 potential (4-body interaction), is explicitly obtained and 
the (statistical) fluctuation is closely explained. The present field the- 
ory model does not conserve energy and is an open-system model. 
(One dimensional) Navier-Stokes equation , Burger's equation, ap- 
pears. In the latter part, we present a way to directly define the dis- 
tribution function by use of the geometry, appearing in the mechanical 
dynamics, and the Feynman's path-integral. 
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1 Introduction 



Boltzmann equation was introduced to explain the second law of the ther- 
modynamics in the dynamical way, in 1872, by Boltzmann. We considers 
the (visco-elastic) fluid matter and examine the dynamical behavior using 
the velocity-field theory. The scale size we consider is far bigger than the 
atomic scale (~ 10 _10 m) and is smaller than or nearly equal to the optical 
microscope scale (~ 10~ 6 m). The equation describes the temporal develop- 
ment of the distribution function /(i,x, v) which shows the probability of a 
fluid-molecule (particle) having the velocity v at the space x and time t. 

We reformulate the Boltzmann equation using the field theory of the ve- 
locity field u(x, '£'). Basically it is based on the minimal energy principle. 
We do not introduce time t. Instead of t, we use the computational step num- 
ber n. The system we consider consists of the huge number of fluid-particles 
(molecules) and the physical quantities, such as energy and entropy, are the 
statistically- averaged ones. It is not obtained by the deterministic way like the 
classical (Newton) mechanics. We introduce the statistical ensemble by using 
the well-established field-theory method, the background-field method® [6]. 
Renormalization phenomenon occurs not because the quantum effect but be- 
cause the statistical fluctuation due to the inevitable uncertainty caused by 
1) step-wise (discrete-time) formulation and 2) continuum formulation for 
the space. 

The dissipative system we consider is characterized by the dissipation 
of energy. Even for the particle classical (Newton) mechanics, the notion of 
energy is obscure when the dissipation occurs. We consider the movement of a 
particle under the influence of the friction force. The emergent heat (energy) 
during the period [ti, t 2 ] can not be written as the following popularly-known 
form. 

/ ^friction dx = ^(^W^WJlt! = E \* ~ E \ti > x i = x ( t i) > x 2 = x {h) , (1) 

where x(t) is the orbit (path) of the particle. It depends on the path (or orbit) 
itself. It cannot be written as the form of difference between some quantity 
at time ti and t 2 . In this situation, we realize the time itself should be re- 
considered when the dissipation occurs. Due to Einstein's idea of "space-time 
democracy" , we have stuck to the standpoint that space and time should be 
treated on the equal footing. We present here the step-wise approach to the 
time-development . 

We do not use time variable. Instead we use the computational-process 
step number n. Hence the increasing of the number n is identified as the time 
development. The connection between step n and step n — 1 is determined 
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by the minimal energy principle. In this sense, time is "emergent" from the 
minimal energy principle. The direction of flow (arrow of time) is built in 
from the beginning. 

In the latter part of this paper, an approach to the statistical averaging 
procedure , based on the geometry of the mechanical dynamics, is presented. 

The content is described as follows. The step-wise dynamical equation 
is presented in Sec. 2. We start with the n-th step energy functional. By 
regarding n steps as the time t n , we derive (1 dim) Navier-Stokes equation. 
In Sec. 3, the orbit (path) of the fluid particle is explained in this step-wise 
formalism. The total energy and the energy rate are also explained. The 
statistical fluctuation is closely explained in Sec. 4. Especially the difference 
from the quantum effect is stressed. Using the path-integral, we take into 
account the fluctuation effect. Owing to the present velocity-field formalism, 
we can obtain Boltzmann's equation, as described in Sec. 5, up to the collision 
term. This step- wise approach is applied to the mechanical system in Sec. 6. 
We take a simple dissipative model: the harmonic oscillator with friction. 
The trajectory is solved in the step-wise way. We find the total energy 
changes as the step proceeds. From the n-step energy expression we can 
extract the geometrical structure (the metric) of the trajectory. The metric 
is used, in Sec. 7, to define the statistical ensemble of the system of N viscous 
particles. We propose some models using the geometrically-basic quantities: 
the length and the area. Conclusion is given in SecJHJ A few appendices 
are provided to supplement the text. App.A treats (1+3) dimensional field 
theory in this step-wise formalism. App.B is the calculation of the statistical 
fluctuation effect and it supplements SecJD An additional mechanical model 
(Spring-Block model) is described with the calculation result in App.C. 



2 Emergent Time and Diffusion (Heat) Equa- 
tion 

We consider 1 dimensional viscous fluid , and the velocity field {u(x); — oo < 
x < oo} describes the velocity distribution in the 1 dim space. Let us take 
the following energy functional [21 [3] of the velocity- field u(x), 

*.[„(*>] -Jm^y + viu) + u ~^~ + l<u - u^f } + n , 

a = l , po = 1 , n = 1,2, • • • , 
1 For a recent review on the nature of time, see ref.pQ. 
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Figure 1: The energy functional I n [u(x)], of the velocity-field u(x). 
In[u(x)] 




' (u (x) } 

Un (x) Un-1 (x) 



v ( u ) = ~l^ u2 + ^y" 4 > it = u{x) , u n —i = it n _i(x) ,(2) 

where po is the constant mass-density and is taken to be 1. H J° is a 'constant' 
term which is independent of u(x). Later we will fix it. a is the viscosity 
constant and is also taken to be 1. El m 2 is the mass density: (the mass of 
the fluid-particle)/ 21. The quantity ([2D is the total energy of the fluid. The 
velocity potential V{u) has the mass term and the 4-body interaction term. 
V 1 {x) is the (ordinary) position-dependent potential. dV d ^ is the external 

2 The time- development term W n {u) is generally written as 

t r / m f , r u ,du. ? dV 1 (x) . T n 

InHx)} = J dx{ — { — ) 2 + V(u) + + W n (u)} + 1° , 

where W„(it) = —(u - u„_i) 2 or - 2w n _i + w n _ 2 ) 2 ■ (3) 

The latter case can be treated in the same way and is given in App.A. Besides the velocity 
field u(x), the fluid matter density field p(x) generally appears. We here consider the 
incompressible case: p(x) = po(constant)= 1. 

We list the physical dimensions of various quantities. [x]=[/]=L, [tt]=L/T, [h]=L/M, 
[m 2 ]=M/L, [A]=MT 2 /L 3 , [cr]=M 2 , [p]= [po]=M/L, [7„]=ML 2 T- 2 , [V}—MLT^ 2 . Further- 
more we have, [%/A^]=(M/L) 3 / 2 T, [Vh 3 \a}= [to~ 3 \/Xa}=[hV>^/m]=T. 

3 We often keep writing a and p~o in order to see the dimensional consistency. 
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source (force) in this velocity-field theory, h is some constant which will be 
identified as the time- separation for one step. u n _i(x) is taken to be a given 
velocity field at the (n-l)-th step. The n-th step velocity field u n (x) is given 
by the minimal principle of the n-th energy functional I n {u). This approach 
is callled "discrete Morse flows method" j2j [3j Hj. 

For simplicity we take the periodic boundary condition for the space. 

u(x) = u(x + 21) , (4) 

where 21 is the periodic length. We may restrict the space region as — I < 
x < I. The variation equation 5I n {u) = 0(u(x) — > u(x) + 5u(x)) gives 

1 a d 2 u n 5V(u n ) dV*(x) 
- Wa; )-« ll _ 1 ( a ;)) = _-— — , 

5V(U) 2 A n 

where we have replaced the minimal solution by u n . From the construction, 
we have the relation: 

4W < 4K-J • (6) 

We, however, cannot say 7 n [w n ] < -^n-i[ w n-i]- The above equation describes 
the n-th step velocity field u n (x) in terms of u n -i(x) and vice versa. Hence 
it can be used for the computer simulation. 

We here introduce the discrete time variable t n as the step number n of 

U n - 

t n = nh = nr x ( — ) , r = h\f\ajm , n — 1, 2, ■ • • 
To 

t = , (7) 

where To is the time unit. The eq.Q is, in terms of the 'renewed' field 
u(x,t), expressed as 

t(u(x, t n _i + h) - u(x, t n _i)) = 
d 2 u(x,t n ) 5V(u(x,t n )) dV 1 (x,t n ) 



dx 2 5u(x,t n ) dx 



4 Note [ro]=T. See the footnote of eq.©. Generally the time t n can be introduced by 
t n = f(n)h where f(n) is a function of n. The form of f(n) defines the time coordinate. 
The change of the form is the transformation of the time coordinate. The simple one is 
f(n) = an + b. In the text f{n) = n is taken. If we take f(n) = —n, the time flow is 
introduced in the inverse way. 
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where we use u(x,t n ) = u n (x), t n = t n _i + h. As h — > 0, we obtain 

du{x,t) _ a d 2 u(x,t) 5V(u(x,t)) _ dV\x,t) 
dt p dx 2 5u(x,t) dx 

where we have replaced both t n and t n -\ by t. This is 1 dim diffusion equation 
with the potential V(u). 

We remind that the variational principle for the n-step energy functional 
I n [u(x)] ©, 5I n = I n [u + 8u] - I n [u] = 0, gives u n (x) for the given u n -i(x). 
We regard the increase of the step number as the time development. H Taking 
into account the fact that, at the (n-l)-step, the matter-particle at the point 
x flows at the speed of u n -i(x), the energy functional I n , (J2J), should be 
replaced by the following one [21 E]. 



I n [u(x)] = J dx {^- o (^) 2 + V(u) + u dV d ^ + ^Hx + hu n -x) - u n ^) 2 } + 7° 

(7=1 , po = l , n = 1,2, • • • 
Tfi^ A 

Note that it(cc) — u n _i(x) in eq.© is replaced by u(x+ hu n _i(x)) — u n _i(x). 
For the simple case of no potential, V = , and no external force, Of- = 0, 

— ) 2 + ^r{ u ( x + ^ "n-l(^)) - Mn-l(x)) 2 } +' Const' ,(11) 

The above functional is equivalent to I n [u(x)} with the potential. 

V(u) = (u(x) - Un ^(x))u n ^(x)^^ + 0(h) . (12) 

where we consider the case of sufficiently-small h. Eq. §5§ gives us the follow- 
ing equation as the minimal equation for J n [u] : 5J n [u] = 

^(u n (x) - u n -i(x)) = - u n -i(x) ^^- + 0(h) . (13) 

Hence the step-wise recursion relation <^ is corrected as 
1 ,du n (x) a d 2 u n 5V(u n ) dV^(x) 

-r(u n (x) -U n ^(x)) +M n _i(x)— = — ■— ,(14) 

h dx po dx 1 bu n dx 



5 Time is defined here by the energy-minimal principle. 
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As done before, let us replace the step number n by the discrete time 
t n = nh. Taking the continuous time limit (h — > 0), we obtain 

du(x,t) du(x,t) a d 2 u(x,t) 5V(u(x,t)) dV l (x,t) 

dt ' dx po dx 2 d~u(x,t) dx 

This is called Burgers's equation (with the velocity potential V(u) and the 
external force ^-) and is considered to be 1 dimensional Navier-Stokes equa- 
tion. Note that the non-linear term in the LHS of eq. ffT5]) appears not from 
the potential ( velocity-field interaction) but from the self-consistency of the 
velocity-field change from step n — 1 to step n. || + w§j = ^ is called 
Lagrange derivative. 

The equation (fl5l) . for the massless case m = 0, is invariant under the 
global Weyl transformation. 

V\x,t) e - 2£ V\e £ x,e 2£ t) , 

u(x,t) — > e~ £ u(e £ x,e 2£ t) , 
d d 

— = d x e e d x , — = d t e~ 2£ d t , 

ox at 

t e 2£ t , x e £ x , (16) 

where e is the real constant parameter. |^| When it is small : |e| 1, the 
variation 5u is given by 

§u = e{-u + x d x u + 2t d t u} + 0(e 2 ) , (17) 

For simplicity, we explain in one space-dimension (dim). The generaliza- 
tion to 2 dim and 3 dim is straightforward. Furthermore the ordinary field 
theory (not using the velocity field but the particle field) is described by this 
step-wise approach in App.A. 



3 Space Orbit (Path) and Total Energy 

The space coordinate x always appears in the velocity field u n (x). We can 
introduce the n-th step space coordinate 

u n {x n ) = , n = 0,1,2,- ■■ 

To 

x n+ i = x n + T u n (x n ) , (18) 

6 If we consider the mass m here appears in some dynamical way (, for example, through 
the spontaneous breakdown ), m — > e~ £ m makes eq. (fT5)l global Weyl invariant. 
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where Xq is a given initial position and r = hy/\a/m in eq.([7j). We can trace 
the position of the matter-point, which was at xq at the initial (0-th) step, 
along the step process: xq, x±, X2, • • ■■ After N steps, the matter-point reaches 

JV-l 

X N = X N _i + T WiV-l(xjV-l) = X + T J2 u n{Xn) ■ (19) 

n=0 

In terms of the continuous time, 

x(T)=x + j-J u(x(t),t)dt , (20) 

where T = N h,t n = n h,x(t n ) = x n ,u(x,t n ) = u n {x). x{t) is the orbit or 
path of the matter-point, at Xq initially, "moving" in the N steps. 
At the n-th step, the total energy of the system, E n , is given by 



f - t r 1 u SIn ^ 
E n = I n [u n \ , where 



Su(x) 



, 



I n [u] = J dx ij~- Q (^) 2 + V(u) + U^- + ^(u(x + hUn-x) - M n _i) 2 } + I° n ,(21) 

The system total energy E n generally changes as the step number increases. 

w[tn) s s 

r dt n 

+ ^{(^+1 + ^—1 «n) - K + hUn-x— 

In ax ax 

where E(t n ) = E n and W(t n ) is the energy rate. From the above formula we 
get the expression for the energy at t = Nh = t n . 

N-l 

E(t N ) = I N [u N ] = Jtv-iK-i] + r G W{t N - X ) = IiK] + r £ W(t n ) . (23) 

n=l 

When we regard the process of the increasing step-number as the time 
development, the system generally does not conserve energy. E(t n ) gener- 
ally changes step by step. We can physically understand that the increase or 



mK+i] - 


- In[u n ] 


TO 




+ (v>n+l 


,dV x 
Un) dx 


-l) 2 }} + 


fo _ fo 

n+1 n ( 2 2) 
T 



7 We will numerically confirm the non-conservation later in Sec|6l 
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decrease of the total system energy is given or subtracted by the outside (en- 
vironment). The energy functional (fTOj) describes the open-system dynamics. 
When E{t n ) satisfies 

W(t n ) = L dE j tn ^ _^ ^(constant) asn^oo , (24) 
To dt n 

we say the system finally reaches the steady energy-state. For the special 
case of Wq = 0, we say the system finally reaches the constant energy-state 

a. 

In SecJHl we treat the Wo = case. As the example of the more general 
case (Wq ^ 0), another model is given in App.C. 

4 Statistical Fluctuation Effect 

We are considering the system of large number of matter-particles, hence the 
physical quantities, such as energy and entropy, are given by some statistical 
average. In the present approach, the system behavior u n (x) = u(x,t n ) 
is completely determined by eq. f[T4"|) when the initial configuration Uq(x) = 
u(x, 0) is given. We have obtained the solution by the continuous variation 
5u(x) to I n [u], ( TTOj) . In this sense, u n (x) is the 'classical path'. Here we 
should note that the present formalism is an effective way to calculate the 
physical properties of this statistical system. Approximation is made in the 
following points: 

1) So far as h ^ 0, the finite time-increment gives uncertainty to the 
minimal solution u n (x). This is because we cannot specify the minimum 
configuration definitely, but can only do it with finite uncertainty. 

2) The real fluid matter is made of many micro particles with small but 
non-zero size. The existence of the characteristic particle size gives 
uncertainty to the minimal solution in this continuum formalism. Fur- 
thermore the particle size is not constant but does distribute in the 
statistical way. The shape of each particle differs. The present contin- 
uum formalism has limitation to describe the real situation accurately. 

3) The system energy generally changes step by step. The present model 
([2]) describes an open-system. It means the present system energetically 
interacts with the outside. Such interaction is caused by the dissipative 
term in ([2]). 

8 Energy constantly comes in or goes out. App.C shows such an example. 

9 Energy does not go out and does not come in. 
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We claim the fluctuation comes not from the quantum effect but from the 
statistics due to the uncertainty which comes from the finite time-separation 
and the spacial distribution of size and shape. 

To take into account this fluctuation effect, we newly define the n-th 
energy functional T[u(x); u n -i(x)} in terms of the original one I n [u(x)}, (fTOl) . 
using the path-integral. 



I n [u(x)} = [ alx{^A 2 + V{u) + u{x) dV + ^r(u(x + hu n ^) - w re _i) 2 } + 1° , 
J 2po ax ax In 

a=l , p = l , V{u) = —u 2 + ^u i ,(25) 

In the above path-integral expression, all paths {u(x); — I < x < 1} are taken 
into account. 

We are considering the minimal path u n (x) as the dominant configuration 
and the small deviation q(x) around it. 



u(x) = u n (x) + y/aq(x) , \y/aq\ < 



SI n [u] 



Su 



, (26) 



In eq.( l2"5"|) and eq.(l2^i). a new expansion parameter a is introduced. 
([a] = [7 n ]=ML 2 T~ 2 ) As the above formula shows, a should be small. The 
concrete form should be chosen depending on problem by problem. It should 
not include Planck constant, h, because the fluctuation does not come from 
the quantum effect. It should be chosen as 

1) the dimension is consistent, 

2) it is proportional to the small scale parameter which characterizes the 

relevant physical phenomena such as the mean-free path of the fluid 
particle, 

3) the precise value should be best-fitted with the experimental data. 

The background-field method [SI E] tells us to do the Taylor-expansion 
around u n . 



/>(*)] = I n [u n (x) + yfaq{x)] = fy/ 2 ^ 

jr^j ll ou[xy 



E S iK\ ,(27) 
z=o 



10 The background-field method was originally introduced to quantize the gravitational 
field theory. Here we borrow the method only to define the statistical distribution measure. 
The present 'splitting' is u = u n + \faq not u = u n + \ffiq. 
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e — r[u n (a;);ii„_i(a:) 



Vqexp 



dx 



q(x) 



Vq(x) exp(—(S + S 2 + 0(q 3 ))) = e -^ Mx)] x 
a 

5 



2! "5u(x) 



) 2 I n [u]\ Un + 0(q 3 ) 
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5u(x 



-In[u]\ Un = 



S = In[u n ] , S 1 = dxq(x) 



8 1 Jul 



8u 



-s 2 

a 



q(x) 2 5 



1 ,dq, 



,m 2 A 



dq. 



1 d , dq. 1 „ ^ , „ a(= 1) d 2 
where we make the Gaussian (quadratic, 1-loop) approximation. 



+ Xu n + m + - 



n-l 



r/.r 



ir[«„(i);u„_i(i)] _ e -±L[u„(x)] x (^ c i e t_D)~ 1 /2 



(detD) 



-1/2 



exp < — Tr In D > = exp < -Tr 



-dr + const 



o r 



(29) 



where r is called Schwinger's proper time[7]. ([r] = [.D _1 ]=L/M.) 

To rigorously define the inside of the above exponent, we introduce Dirac's 
abstract state vector \x > and < x\. 



D 



a(=l) d 2 
p (= 1) dx' 



+ D)G(x,y-r) = 
- V(x) , V(x) -- 



<x\e TD \y>=G{x,y,T) 
]xm o G(x,y;r) = 5(x - y) 

2 __„2 1 dll n —l 



-Xu n (x) - m - - 
n 



dx 



,(30) 



G(x, ?/; r) is called heat-kernel. 

In App.B, we evaluate ln(det .D) 



'1/2 _ I 



SS°dVTrG(x,y)/T 



\ Jq 00 dr /I; <ix G(x, x)/r. Up to the first order of V, the result is given by 



I ri 1 
= VeA - - / dze(Xu n (z) 2 + m 2 + - 

it 2y/7iefi J-i h 



dz 



(31) 



where the infrared cut-off parameter u = \/(x/Z and the ultraviolet cut-off 
parameter A = h~ l are introduced, o e _1 = a/p = 1. We see the mass 

11 0(/j) = (h/2)iin_i(dq/dx) 2 may be ignored for ft<l. 

12 The last term in the round-brackets does not contribute due to the periodicity of 
u n -x(z). 
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parameter m 2 shifts under the influence of the fluctuation. 



m 2 ->■ m 2 + — =e\ = m 2 + aAW ^= , (32) 
And the bottom of the potential shifts as 

) -»■ V(M m 




+ m /-p- {-=- J— (m 2 + -)} . (33) 



The coupling A is also shifted by the 0(V 2 ) correction. The shift of 
these parameters corresponds to the renormalization in the field theory [9]. In 
this effective approach, we have physical cut-offs \x and A which are expressed 
by the (finite) parameters appearing in the starting energy-functional. When 
the functional (I1UI) (effectively) works well, all effects of the statistical fluctu- 
ation reduces to the simple shift of the original parameters. This corresponds 
to the renormalizability condition in the field theory. 



5 Boltzmann's Transport Equation 

We use, for simplicity, the the original names for the shifted parameters. The 
step-wise development equation ( TH|) with SV/Su = m 2 u + ^u 3 + M n _i^-, 
= 0, and the correction term (fT2l is written as 

1 / w d o ^3 dll n 

-r{Un{x) - it n _i(x)) = — -T- - m u n - — u n -M n _i— — or 
h ax 1 3! ax 

dx 



13 This corresponds to renormalization of "mass" m 2 in the field theory. When natural 
cut-offs (IR and UV) are there in the system model-parameters, the divergences coming 
from the space integral and the mode summation are effectively expressed by "large" but 
finite quantities. 

14 The coupling (A) shift can be obtained from 0(V) result ([55)) by assuming the "renor- 
malization" consistently works. Noting 7(u mm ) = — 6m 4 / A, A should shift as 



a t A , 2 / Po 1 J , I I 



6 m? V 47T(t y/h I y \foh 
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The latter form is convenient for the 'backward' recursive computation: u n — > 
u n _i. When the system reaches the equilibrium state after sufficient recursive 
computation (n 3> 1), we may assume u n -\(x) = u n (x) = u°°(x). u°°(x) 
satisfies 

d 2 n°° \ rlii 00 

%-^°°-^-^ = • < 35 > 

We here introduce the distribution function f n (x, v) as the probability for 
the matter-point particle in the space interval x ~ x + dx and the velocity 
interval v ~ v + dv, at the step n, is given by 

-=-f n (x,v)dxdv , (36) 

where N n is the total particle number of the system at the step n.[8] Then 
the n-th distribution f n (x,v) and the equilibrium distribution f°°(x,v) are 
introduced as 

U°°(x) = -r^ / Vf°°(x,v)dv , U n (x) = — -r-- [ Vf n (x,v)dv , 

u n (x) — > u°°(x) and / n (x, v) — > f°°(x, v ) as n — >• oo , (37) 

where m°°(x) is the equilibrium velocity distribution. p n {x) is the particle 
number density. The continuity equation is given by 

\{pn{%) ~ Pn-l(x)) + ^-{p n {x)u n {x)) = . (38) 

n dx 



The recursion relation (1341) is expressed, in terms of the distribution func- 
tions, as 

r + hu n -l{x), V) - f n -i(x, v)] = 

9 2 fn(x, v) 2 A . 2 
— m f n [X, V) - gj Jn(2G, U)«n(«) , 

where m„(x) = — ^— r ^ vf n (x,v)dv , (39) 

This is the Boltzmann's transport equation for the 2-body and 4-body velocity- 
interactions. We can express the step-wise expression (1311 in the continuous 



15 Here we list the physical dimension of some quantities appearing this section. [a;]=L, 
K] = H=L/T, [fc B 7l=ML 2 /T 2 , [Pn(x)]=L-\ [p^x^ML- 1 , [/„]=T/L 2 , [P„]=ML/T 2 , 



[g n ]=ML 2 /T 3 
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time t form as in Sec|2J This is the integrodifferential equation for f n (x,v) 
when p n is a constant. The right hand side (RHS) is called collision term. 
We notice when we may replace w n _i, in the LHS of eq. (|39l . by u n , the above 
recursion relation determine the (n-l)th step distribution f n -i by the n-th 
step data, /„ and u n . El 

We now introduce some physical quantities used in the non- equilibrium 
statistical mechanics. The entropy S n and the total particle-number N n are 
defined by 

S n = ~ks J dv J dxf n (x,v)\nf n (x,v) , 
N n = dxp n (x) = dx dvf n (x,v) , p n (x) = / dvf n (x,v) , (40) 



where ks is Boltzmann's constant. The momentum conservation at each 
point, x, requires 

= p n (x) f dv(v - u n (x))f n (x, v) , u n (x) = f dvv f n (x, v) .(41) 

Besides the particle- number density p n (x), we have introduced the mass den- 
sity p n (x) at step n. 

We here consider the case of one kind particle. 

^4 = m i ( constant ) , m x = m 2 x (2/) , (42) 

where m\ is the particle mass. El In this case the total mass M n is given by 

M n = J dxp n (x) = m 1 J dxp n (x) = miN n . (43) 

The temperature distribution T n (x), the heat current distribution q n (x) and 
the pressure P n (x) are given by 

\k B T n {x) = ■ -r^ f dv^(v - U n (x)) 2 f n (x,v) , 

2 pJx) J 2 



Pn(x) 

q n (x) = J dv^-(v - u n {x)ff n {x,v) 
P n (x)=m 1 dv(v-u n (x)) 2 f n (x,v) , (44) 



16 Solving the recursion equation in the backward way, with respect to the step 
number n, is convenient in practice. 

17 In Sec[2j we considered the case of a = 1 and p = po = 1 for simplicity. When we 
generalize it, we have to use some equations appearing in this section to determine the 
form of a n {x) and p n {x). Especially the viscosity a n {x) depends on the temperature T n {x) 
introduced in (|33|). 

1 For k kinds particles, we introduce p n and p n fields for each one: {p n l , ph\i = 
1, 2, • • • , k} with the relation p l n [x) / p n l (x) = 
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Figure 2: The harmonic oscillator with friction, (j49p . 
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where ks is Boltzmann's constant. 

Using the transport equation which f n (x,v) satisfies, we expect the 
following two equations are satisfied. The first one is the rephrasing the 
Navier-Stokes equation ([3] 



^(U n (x) - tt„_i(x)) + M n-1^7 

The second one is the energy equation. 



1 dP n (x) 
p n (x) dx 



(45) 



\p n {x)k B {hr n {x) - r n - x {x)) + u n J^} = _ p n (x)^^- .(46) 

2 n dx dx dx 

In the remaining sections, we present an alternative approach to the dis- 
tribution function f n (x,v). 



6 Classical and Quantum Mechanics and Its 
Trajectory Geometry 

We can treat the classical mechanics and its quantization ( the quantum 
mechanics, not the quantum field theory) in the same way. In this case, the 
model is simpler than the previous case (space-field theory) and we can see 
the geometrical structure clearly. Let us begin with the energy function of 
a system variable , x, (1 degree of freedom). For example the position (in 
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1 dimensional space) of the harmonic oscillator with friction. We take the 
following n-th energy function to define the step flow. 



71 Tfl 

K n (x) = V(x) + 1 ^(x-x n ^ 1 ) 2 + —(x-2x n „ 1 +x n _ 2 ) 2 + K° , (47) 

where V(x) is the general potential and K® is a constant which does not 
depend on x. For the harmonic oscillator V(x) = kx 2 /2 where k is the 

§ring constant, r] is the friction coefficient and m is the particle mass. 
We assume x n _i and x n _ 2 are given values. As in Secf2J the n-th step 
position x n is given by the minimal principle of the n-th energy function 
K n (x): SK n = 0, x — > x + Sx. 



sv 

Sx 



Tj Tfl 

+ - x„_i) + — (x n - 2x n _i + x„_ 2 ) = , (4E 
_ n n z 

X — Xn 



With the time t n ([7J, the continuous limit (h — > 0) gives us 

dV(x) dx d 2 x , . 

where t n = nh — >• t, x n = x{t n ) x(t), (x n — x n _i)/h = dx/dt\ tn — >■ 
dx/dt, (x n — 2x„_i + x n -2)/h = d 2 x/dt 2 \ tn —> d 2 x/dt 2 . For the case of 
V = kx 2 /2, this is the harmonic oscillator with the friction r]. See FigfJl 
This is a simple dissipative system. 

The recursion relation fj4*8|) gives us, for the initial data Xq and x±, the 
series { x n = x{t n )\n = 0, 1,2, ■ • • }. This is the classical 'path'. The fluc- 
tuation of the path comes from the uncertainty principle of the quantum 
mechanics in this case. (We are treating the system of 1 degree of freedom. 
No statistical procedure is necessary. ) As the time-interval h tends to zero, 
the energy uncertainty grows (At ■ AE > h). Hence the path x n , obtained by 
the recursion relation (1481) . has more uncertainty as h goes to 0. As in SecJH 
we can generalize the n-th energy function K n (x), (1471) . to the following one 



19 Here we list the dimension of parameters and variables in this section. [x] = [x„]=L. 
H = K]=L/T, [t]=[t n ]=T, [ (7 ]=T 1 /2 M - 1 /2 j [?i] = ML 2 /T, [m]=M, [r?]=M/T, [h]=T, 
[if„] = [y]=ML 2 /T 2 , [/c]=M/T 2 , [V^^M 1 / 2 , [%/m/i 2 ]=M 1 / 2 T. Some ones, such as 
t, m, h, q and V, appearing before this section have different dimensions in this section. 

20 The eq. (T4^| is compared with the eq. © for the case of the no external force: 

_ SV(u(x,t)) ^ d 2 u(x,t) . du(x,t) () 
Su(x,t) dx 2 ^° dt ' 

where we notice the friction term in ea.(|49p corresponds to the dissipative term in eq.©. 
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r(x n _i,x n _ 2 ) in order to take into account the quantum effect. 

e -^hr(xn-i,xn_ a ) _ f°° e~T\ hKn< y x ) 
J — oo 

Ti Tfl 

K n (x) = V(x) + ^{x-x n - 1 ) 2 + —{x-2x n - 1 +x n - 2 ) 2 + K n ' . (50) 

We can evaluate the quantum effect by the expansion around the classical 
value x n : x = x n + Vh q where % is Planck constant. El 



e -±hT(x n ;xn-i,x n - 2 ) _ J ^ x e ~\hK n (x) _ j ^ e -i/i_ftT n (a: n +%) 



5x 



5]/ 

5x 



?7 m 

3? — n. 



Tl Tj Tfl 

Y n = T(x n ; x n _i, a; n _ 2 ) = A' n (x n ) + — ln(/c + - + — ) , (51) 

where the final expression is for the oscillator model: V = kx 2 /2. The 
quantum effect does not depend on the step number n. It means the quantum 
effect contributes to the energy as an additional fixed constant at each step. 
The energy rate is obtained as 

dK(t n ) _ ( \ jt- t \ -p r _ , dr(t n ) 

a T. = l\ n +l{Xn+l) — J^n\X n ) — 1 n+1 — 1 n = I- 



= V(x n+ i) - V(x n ) + ^{(l n+ i - X n ) 2 - (X n - X n -i) 2 } 

Tfl 

+ ^{(x n+1 -2x n + x n - 1 ) 2 - (x n -2x n - 1 + x n - 2 ) 2 }K° n+1 - K n . (52) 

The present system is again an open system, and the energy generally changes. 

In terms of the position difference x n — x n _\ = Ax n and the velocity 
difference (x n — 2x n „ 1 + i n _ 2 )//i = v n — v n _i = Av n , we can rewrite the 
energy at n-step and read the metric as follows. El 

Tj Tfl 

K n (x n ) = V(x n ) + — (i„ - x n _i) 2 + — (a; n - 2x n _i + x n _ 2 ) 2 + 

= l 2 {V(x n )(At) 2 + f(Ax n ) 2 + ^(Av n ) 2 } + K° n , (53) 

where h (time increment) in the first term within the round brackets is re- 
placed by At. This shows us the metric for the n-step energy function is 
given by 



ivL) 2 , 



(As n ) 2 = 2h 2 K n (x n ) = 2V(x n '/^h)(At) 2 + (Ax n ') 2 + (At 

x n ' = \/r]hx n , vj = V mh 2 v n , (54) 



21 Do not confuse h (Planck constant/27r) with h (time interval). 

22 v n = (x n - x n -x)/h = Ax n /h 
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where , for the oscillator model, V(x n '/ y/rjh) = (k'/2)x n ' 2 , k! = kjr\h. Eg. ( I54|) 
shows the energy-line element As 2 in the (t, x n \ v n ') space. El Note that the 
above metric is along the path x n = x(t n ), v n = v{t n ) = (x(t n ) — x(t n _i))/h 
given by (HHl) . The metric is used, in the next section, as the geometrical 
basis for fixing the statistical ensemble. 

We take the freedom of the value K® in the following way. 



Kl = -V{x n ) - — {x n - 2x n _! + x„_ 2 ) 2 + V{x ) + —3(2:1 - x ) 2 . (55) 



This is chosen in such a way that the step n energy K n (x n ), for the no 
dissipation (77 = 0), does not depend on the step number n and the value 
is the total energy at the initial step (last 2 terms in (155]) ). The graphs 
of movement and energy change, for various viscosities, are shown in FigJS]- 
[HJ For the no friction case, the oscillator keeps the initial energy (Figfjj). 
When the viscous effect appears, the energy changes step by step, and finally 
reaches a constant nonzero value (Figj6jFigJ71 Figj9]). We understand the 
finally-remaining energy (constant) as the dissipative one. Physically (in 
the real matter) the energy is realized as the pressure and the temperature 



which characterize the particle's "environment" (out-side world). For the 
resonate case (4k/m = (77/m) 2 ), both the movement and the energy are 



7 Statistical Ensemble, Geometry and Initial 
Condition 

In this section, we consider some statistical ensemble of the classical me- 
chanical system taken in the previous section. Namely, we take iV 'copies' 
of the classical model and regard them as a set of (1 dimensional) particles, 
where the dynamical configuration distributes in the probabilistic way. N is 
a large number. The set has N degrees of freedom: xi, x 2 , ■ ■ ■ , x^. As the 
physical systems, (1 dimensional) viscous gas and viscous liquid are exam- 
pies. Each particle obeys the (step -wise) Newton's law (I48j) with different 
initial conditions. N is so large that we do not or can not observe the ini- 
tial data. Usually we do not have interest in the trajectory of every particle 

23 In ca. flMl) . the first term shows the potential part, the second one the kinetic part 
and the third one a new term. In ref.|10) and ref.|ll). the hysteresis term appears as a 
new one. 

24 Note that the particle does not move at the final stage. 

25 For example, N - 10 23 . 

26 We are considering ./V-body problem where each particle moves (fluid flows) with 
friction. We aproach it using the effective 1-body energy function (H71) . 




large. 
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Figure 3: Harmonic oscillator with no friction, Movement, 

h— 0.1, ^Jk/m— €.01, rj/m— 0.0, total step no =20000. The step-wise solu- 
tion (jJSJ) correctly reproduces the analytic solution: x(t) = sin(o;t) , < 
t < 2000 , x(0) = 0, x(0) = 1. 
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Figure 4: Harmonic oscillator with no friction, Energy change, 

h=0.1,yjk/m=0.01, r]/m=0.0, total step no =20000. 
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Figure 5: Harmonic oscillator with friction, Movement, /i=0.1,y /c/m=0.01, 
77/771=0. 005 (Elasticity dominate) and 0.03 (Viscosity dominate), total step 
no =20000. The step- wise solution (148]) correctly reproduces the analytic 
solution: (l)Elast.Dom. x(t) = e -V*/2 s i n (y / 4o; 2 - t/ 2 t/2) , < t < 
2000 , x(0) = 0, x(0) = ^Au 2 - r]' 2 /2 = (1.94 x 10- 2 )/2 (2)Visco.Dom. 

< t < 2000 , x(0) = 0, x(0) = 



x (t) = e -V*/2 S inh(^r/' 2 -4w 2 t/2) 
4w 2 /2 = (2.24 x 10~ 2 )/2. 
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Figure 6: Harmonic oscillator with friction, Energy change, 

h=0.1,^k/m=0.01, r?/m=0.005 (Elasticity dominate) , total step no =20000. 
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Figure 7: Harmonic oscillator with friction, Energy change, 

h— 0.1, yjk/m— 0.01, i]/m=0.03 (Viscosity dominate), total step no =20000. 



Q_ 
CD 

-I—' 

C/) 
i 

C 



D) 
i_ 
CD 
C 
LU 



c 
X 



0.00012474 



0.00012472 



0.0001247 



0.00012468 



0.00012466 



0.00012464 



0.00012462 



0.0001246 



"EnefrHO(2).dat" 



Viscosity Dominate 



h=0.1, om=0.01, eta/m=0.03 



5 10 15 

n(Step No) x h(Step Intvl) x om(Freq) 



20 



23 



Figure 8: Harmonic oscillator with friction, Movement, h=0.1,yk/m=0.01, 
rj/m=0.02 (Resonant), total step no =20000. The step-wise solution ( |4"8"]1 
correctly reproduces the analytic solution: x{t) = te _aJ * , < t < 
2000 , x(0) = 0, x(0) = 1. 
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Figure 9: Harmonic oscillator with friction, Energy change, 

h=0.1,^k/m=0.01, r]/m=0.02 (Resonant), total step no =20000. 
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and do not observe it. We have interest only in the macroscopic quantities 
such as total energy and total entropy. The N particles (fluid molecules) in 
the present system are weakly interacting each other in such way that each 
particle almost independently moves except that energy is exchanged. 

As the statistical ensemble, we adopt the Feynman's idea of "path-integral" 
H21 H31 [HI HS1 [13 H3 HH |in|. We take into account all possible paths {y n }. 
{y n } need not satisfy (118]) nor certain initial condition. As the measure for 
the summation (integral) over all paths, we propose the following ones based 
on the geometry of ( |5l|) . As the first measure, we construct it in terms of 
the length, using the "Dirac-type" metric [TT| [TO]. 

(ds 2 ) D = 2V(X)dt 2 + dX 2 + dP 2 - on-path(X = y(t), P = w(t)) -> 

(2V(y) + y 2 + w 2 )dt 2 , 

r? rP i p/h i 

L D = ds\ on _ path = <j2V(y) + y* + w*dt = h^y/2V(y n ) + yl + w* , 
Jo Jo n=0 

dfi = e~i LD VyVw , e _/3F = f \{dy n dw n e~^ LD ,(56) 

where a is a parameter with the dimension of length ([a]=L). See FigJTQl 
As explained in SecJH it is appropriately chosen problem by problem. (3 is 
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introduced to restrict the t-region (0 < t < fi) and, in this context, should 
be regarded as a part of the choice of the ensemble. /3 plays the role of the 
inverse temperature. Among all possible paths {y n }, the minimal length 
(5Ld = 0) solution, (H8|) . gives the dominant path {x n }. 

The second choice is constructed using the "Standard-type" metric. 

(ds 2 ) s = ^[(ds 2 ) D ? ~ on-path(X = y(t),P = w(t)) - 



(2V(y) +f + w 2 fdt 2 

L S = J ds\„ ath = J (2V{y) + y 2 + w 2 )dt = h^V{y n )+y 2 n + w 2 n ) 



n=0 
0/h 



>±(2V(y n )+yl) 



dfi = e <* Ls VyVw , e l3F = / \\_dy n dw n e <* Ls = (const) / dy 

where we should notice dt (= h > 0) is non-zero. In both cases above we 
take the metric of the 3 dimensional (bulk) space-time (X, P, t) , which is 
consistently chosen with the trajectory metric (1541) . Note that the standard 
case is exactly the same expression as the free energy (trace of the density 
matrix) expression in the Feynman's textbook[18j. 

Another choice of path is making use of surfaces instead of lines. Let us 
consider the following 2 dim surface in the 3 dim manifold (X, P, t). 

X 2 + P 2 = r 2 (t) , 0<t<(3 , (58) 

where r(t) is arbitrary (non-negative) function of t. We respect here the 
isotropy of the 2 dim phase space (X, P). See FigJTTl By varying the form of 
{r(t) : < t < f3}, we obtain different surfaces. Regarding each of them as a 
path used in the Feynman's path-integral, we obtain the following statistical 
ensemble. First the induced metric gij on the surface (|58i) is given as 

^ 2 )-lon-path= 2 ^W^ + rfX2 

(9ij 

where (X X ,X 2 ) = (X,P). Then the area of the surface 



A = J ^det g tj d 2 X = J J 1 + ^jdXdP 




27 fi/h — N should be an (large) integer. The increment h is the (inverse) temperature 
unit as well as the time unit. From the dimensional analysis kB l rjl 2 /l3 corresponds to 
the temperature. Here is Boltzmann's constant and rjl 2 = ft is the combination of the 
friction coefficient r\ and some length scale I (|7]=L) such as the mean free path of the fluid 
particle. Note that fi has the same dimension as Ti. [?i'] = [?i]=ML 2 /T. 
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Figure 11: The two dimensional surface, ( 155]) . in 3D bulk space (X,P,t). 

(X=x (t) , P=p (t) ) . ' 




We consider all possible surfaces of ( 155]) . The statistical distribution is, using 
the area A, given by 

= J™ dp J r{Q) = p Y[VX(t)VP(t)e^ A , (61) 

r(f3) = p 

In relation to Boltzmann's equation (SecJ5J, we have directly defined the 
distribution function f(t,x,v) using the geometrical quantities in the 3 dim 
bulk space. 



8 Conclusion 

We have presented the field theory approach to Boltzmann's transport equa- 
tion. The collision term is explicitly obtained. Time is not used, instead the 
step number n plays the role. We have presented the n-th energy functional 
(ITU]) which gives the step n configuration u n (x) from the minimal energy 
principle. We regard the step flow ( the increase of n ) as the evolution of 
the system, namely, time- development. Navier-Stokes equation is obtained 
by identifying time t as nh (|7]). Fluctuation effect due to the micro struc- 
ture and micro (step-wise) movement is taken into account by generalizing 
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the n-th energy functional I n [u(x)], fTTU]) . to T[u(x); u n _i(x)}, ( 12"oT) . where the 
classical path u n (x) is dominant but all possible paths are taken into account 
(path- integral) . Renormalization is explicitly done. The total energy gener- 
ally does not conserve. The system is an open one, namely, the energy comes 
in from or go out to the outside. In the latter part we have presented a 
direct approach to the distribution function f n (x,v) based on the geometry 
emerging from the mechanical (particle-orbit) dynamics. We have examined 
the dissipative system using the minimal (variational) principle which is the 
key principle in the standard field theory [T9]. 



9 Appendix A (3+l)D Scalar Field Theory 

3+1 dimensional scalar field is here treated in the present step- wise formal- 
ism. We start with the following n-th step energy functional. 

/ n [0(x)] = J a?x{ l -{V<P) 2 + V((f>) + ^(0 - 20 n „! + n _ 2 ) 2 } n = 2, 3, • ■ • 



m , 2 A (4 



v{4>) = -j(j> + ^ 



(x) 



= <pn-l(X), n _ 2 



2 X 



where n ,_2 and 4> n -i are given, (x) = (xi,x 2 , x 3 ) is the 3 dimensional spacial 
coordinates. The step n configuration n (x) is defined to be the energy 
minimal one. 



6L 



n2j SV 



+ 7^(0n - 20 n _i + n _ 2 ) = 

^ = m0+ 3! 



(63) 



Using the step-time notation: n (x) = 0(x, t n ), t n = nh, we obtain, in the 
continuous time limit h — > +0, 



- V 2 + m 2 )0 + A 3 = o 



(64) 



This is the (3+1) dim 4 scalar field equation. 



10 Appendix B Calculation of Fluctuation 
Effect 

In SecJH we have developed the method of calculating the statistical fluc- 
tuation occurring in the (1-dim) viscous fluid matter. The background-field 
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method is taken. At the 1-loop approximation, the key quantity to calcu- 
late the energy functional is the heat-kernel G(x,y;r) given by (see 
eq.fQD) 

< x\e~ TD \y >= G(x,y;r) , 

d 

(— + D)G(x,y;r) = , lirn G(x, y; r) = 6(x - y) , 
or r^+o 

D = -e-^-V(x) , V{x) = -\u n {xf- m *-\- d -^ , (65) 

where \x > and < x\ are Dirac's abstract state vectors, e = p^jo is explicitly 
written to show the dimension consistency. In the text, we take e = 1. For 
the calculation, in this appendix, we change the scale of r and D as follows. 

r -> - = f , D ->■ eD = D , 
e 

D = -^- 2 -eV{x) , [f]=L 2 , [,D]=L- 2 , [e]=l/LM . (66) 

In the following within this appendix, for simplicity we omit the symbol of 
'tilde'. 

The kernel is formally solved as 
G(x, y; t) = Gq(x — y; r) + J dz J duS(x — z;t — u)eV(z)G(z,y;u) ,(67) 
where Go(x — y; r) and the (heat-) propagator S(x — y; r) are given by 

ff.( I -^)=r^e- ft+ " | -" = ^- 6 # , r>0 , 
j-oo An V47rr 

, f f°° dk dk° exp{-ik°r + zfc(:r - y)} 

They satisfy 
d d 2 

(-Q--g-^) s (x-y',V^ s ( T ) S ( x -y) > ^^(x -y;r) = 5(x-y) .(69) 
Up to the first order of V", 

G(x, x] t) = -== + fdz f T duS{x -z;t- u)eV{z)G {z -x;u) + 0{V 2 ) .(70) 
V Atit J Jo 
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Figure 12: The spring-block model, (I75p . 
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The second term is evaluated as 

J dz J duGo(x — z;t — u)eV(z)Go(z — x; u) = 

fdzTdu^- , 1 eV(z)exp{- . T (x - zf} . (71) 



Finally the contribution to ln(det D)- 1 ' 2 = \ / oo (dr/r)TrG(x, y; r) = \ / °°(dr/r) ti dxG(x, x; r) 
is evaluated as 

2 1 dT 2V^I dzf9iz) = ^l -- 1/2 */^W = 

where the infrared cut-off parameter ji = -y/o"// an d t ne ultraviolet cut-off 
parameter A = h~ l are introduced. 



11 Appendix C Spring-Block Model 

In Secj6j the movement of the harmonic oscillator with friction was examined. 
Here we treat the movement of a block which is pulled by the spring which 

28 The dimensions of these parameters are [/t] = [A]=M/L. The space- integral part 
(/ dx-- •) in (|72| is evaluated as /_ ; dx exp{— (t/4(t — oj)u)){x — z) 2 } ~ J^expj"} = 
2 \J 7r(r — lj)l>j/t where I is safely extended to infinity. 
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moves at the constant speed V. The block moves on the surface with friction. 
We take the following n-th energy function to define the step flow. 

— 77 77) 

K n {x) = V(x) - hnkVx + -Ux - x„_!) 2 + —{x - 2x n _ 1 + x n _ 2 ) 2 + K° n , 

V(x) = ^ + klx ,(73) 

where t] is the friction coefficient and m is the block mass. The potential 
V(x) has two terms: one is the harmonic oscillator with the spring constant 
k, and the other is the linear term of x with a new parameter / (the natural 
length of the spring). V is the velocity (constant) with which the front-end 
of the spring moves. K® is a constant which does not depend on x. It will be 
fixed later. The n-th step x n is determined by the energy minimum principle: 
SK n (x)\ x=Xn = 0. 

— (x n + I - nhV) + H — -(x n - x n _i) + -rx(x n - 2x n _i + x n _ 2 ) =0 or 
m m n hr 

_ u 2 {-l + nhV) + ^x w _i + ^2 (2x n _i - x n _ 2 ) _ jk , _ rj 
u ; + 1 + F » m m 

For the continuous limit: /i — > 0, n/i = t n — > t, (x n — x n -\)/h — > x, (x n — 
2x n -i + x n -2)/h 2 — > x, the above recursion relation reduces to the following 
differential equation, o 

mi = k{Vt — x — I) — rjx . (75) 

We keep the step- wise approach. The system energy given by K n (x n ). 
Taking the constant term K® as 

777 777 
Kl = -V(x n ) - 7^2-(Xn - 2X„_! + X n „ 2 ) 2 + V(X ) + — (s a - X ) 2 

-t-Zm/cl^Xn + hkVx n , (76) 

the energy is given as 

K n (x n ) = hkVx n + ^(x n - x^) 2 + V(x ) + ^(xx - x ) 2 , (77) 

We have taken the constant term K®, (176 p . in such a way that the system 
keeps the constant energy when the energy dissipation does not occur (r] = 
0, V — 0). The last two terms in ( ITS]) comes from the following relation. 

(^kx 2 - Wxt + = -7]i 2 - kVx . (78) 



29 This equation is called spring-block model and is used to explain some aspect (stick- 
slip motion, etc) of the earthquake. 
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Note that the y-terms appear both hand-sides, the total derivative terms 
and the dissipative terms. I 30 l 

The graphs of movement (x n , eq. (1T4"l) ) and energy change (K n (x n ), eq. (f77|) ) 
are shown in Fig{13] and FigffJ] respectively. From the graph of FigJT4"l we 
see this system reaches the steady energy-state as n — > oo, ( ]24"j) . 
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Figure 14: Spring-Block Model, Energy Change, /i=0.0001,y / /c/m=10.0, 
77/m=1.0, V=1.0, T=1.0, total step no =20000. 
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